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The distribution of gluon fields in hadrons is of fundamental interest in QCD. Using 
lattice QCD we have observed the formation of gluon flux tubes within tri-quark (baryon) 
systems for a wide variety of spatial distributions of the color sources. In particular we 
have investigated configurations where two of the quarks are close together and the third 
quark is some distance away, which approximates a quark plus diquark string. We find 
that the string tension of the quark - diquark string is the same as that of the quark - 
antiquark string on the same lattice. We also compare the longitudinal and transverse 
profiles of the gluon flux tubes for both sets of strings, and find them to be of similar 
radii and to have similar vacuum suppression. 



1. Introduction 

Recently there has been renewed interest in studying the distribution of quark and 
gluon fields in the three-quark static-baryon system. While the earliest studies were 
inconclusive^, improved computing resources and analys is tech niques now make it 
possible to study this system in a quantitative manne J^j^t In particular, it is 
possible to directly compute the gluon field distributio n 1 ^ ! ^ using lattice Q CD 
techniques similar to those pioneered in mesonic static-quark systemd ^O I H I We 
first investigate the formation of flux tubes, where vacuum gluon field fluctuations 
are suppressed, in systems where the three quarks are approximately equidistant. 
In order to compare baryonic and mesonic flux tubes, we have further investigated 
three-quark static baryons where two of the quarks are close together and the third 
is some distance away 
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Fig. 1. Gauge-link paths or "staples," Ui, U2 and U3, forming a three-quark Wilson loop with the 
quarks located at r\, r*2 and r*3. e abc and e a b c denote colour anti-symmetrisation at the source 
and sink respectively, while r indicates evolution of the three-quark system in Euclidean time. 



2. Baryon Flux tubes on the Lattice 

In order to study flux-tubes on the lattice, we begin with the standard approach 
of connecting static quark propagators with spatial-link paths in a gauge invariant 
manner. We use APE-smeared spatial-link paths (with a — 0.7 and 30 smearing 
steps in all cases) to propagate the quarks from a common origin to their spatial 
positions as illustrated in Fig. Q] The static quark propagators are constructed 
from time directed link products at fixed spatial coordinate, YliU t {x,ti), using 
the untouched 'thin" links of the gauge configuration. Finally smeared-link spatial 
paths propagate the quarks back to the common spatial origin. The three-quark 
Wilson loop is thus defined as: 

W 3Q = l £ abc £ a ' b ' c ' Jjf Uf Uf ', (1) 

where Uj is a staple made of path-ordered link variables as shown in Fig. [1] We 
consider Y and T shaped paths in the x — y plane, where the three quarks are created 
at the origin, then are propagated to positions r*i, r 2 and r 3 , which approximate 
an equilateral triangle, before being propagated through time and finally back to a 
sink at the same spatial location as the source. We investigate seven such triangles, 
with average inter-quark distance varying from 0.27 fm to 1.7 fm. 

The gluon-field fluctuations are characterised by the gauge-invariant action den- 
sity S(y,t). We calculate the action density using the highly-improved 0(a 4 ) three- 
loop improved lattice field-strength tensor on four-sweep APE-smeared gauge 
links. 

For a Wilson loop of Euclidean time extent r we evaluate the correlation function 

C{y;r 1 ,r 2 ,r 3 ;T) = . - - — . (2) 

{W 3 Q{r 1 ,r2,r 3 ;T)) {S{y,r/2)) 

where (• ■ • } denotes averaging over configurations and lattice symmetries. This cor- 
relates the quark positions with the gauge-field action in a gauge invariant manner, 
and has the advantage of being positive definite, eliminating any sign ambiguity on 
whether vacuum field fluctuations are enhanced or suppressed in the presence of 
static quarks. For y well away from the quark positions, there are no correlations 
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Fig. 2. Expulsion of gluon-field fluctuations from the region of static quark sources illustrated by 
the spheres. An isosurface of C(y) is illustrated by the translucent surfaces. The surface plots (or 
rubber sheets) describes the values of C(y) for y in the quark plane, (y\, y2,0). Results are for 
30-sweep smeared T-shape (left) and Y-shape (right) sources with the largest quark separation 
considered. 



and C — *■ 1. We find that C is generally less than 1, signaling the expulsion of 
vacuum fluctuations from the interior of heavy-quark hadrons. 

We first consider 200 quenched QCD gauge-field configurations created with the 
C(a 2 )-mean-field improved Luscher-Weisz plaquette plus rectangle gauge actiorP31 
on 16 3 x 32 lattices at (3 = 4.60. The long dimension is taken as being the x 
direction making the spatial volume 16 2 x 32. Using a physical string tension of 
er = (0.440 GeV) 2 = 0.981 GeV/fm, the QQ potential sets the lattice spacing to 
a = 0.123(2)fm. 

We use lattice symmetries to improve the signal to noise ratio of our measure- 
ments. These include translational invariance (any point on the lattice can be taken 
as the origin), reflection in the x plane and 90° rotational symmetry about the x- 
axis. The advantage of this approach is that we do not have to perform any gauge 
fixing to find a signal in the flux distributions. 

In Fig. [2] we show examples of the expulsion of vacuum fluctuations and the 
formation of flux-tubes for our largest T and Y shaped configurations. 

A detailed analysis of these flux tubes is given elsewher e!^ !. The major con- 
clusions can be summarised. We do not see any evidence for A - shaped (empty 
triangle) flux tubes. At large quark separations we do see Y shaped flux tubes, even 
from initial T shapes, which relax towards a Y shape. We observe a potential which 
is a linear function of the minimum length of string needed to connect the quarks 
to the Fermat point, and the extracted string tension (a = 0.97(1) GeV fm _1 ) is 
in good agreement with the quark - antiquark string tension. 



3. Diquarks on the Lattice 

There has recently been a renewal of interest in the properties of diquarks in 
hadronic systems, as they may play an important role in the existence of ex otic 
states, such as the putative 9 + , or in explaining the scarcity of such exoticJ^. 
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Fig. 3. Expulsion of gluon- field fluctuations from the region of static quark sources illustrated 
by the spheres. An isosurface of C(y) is illustrated by the translucent surface. A surface plot (or 
rubber sheet) describes the values of C(y) for y in the quark plane, (y\, j/2i 0). 

In QCD, two quarks close together, a diquark, can transform either according to 
the conjugate representation (3) or the sextet (6) representation of SU(3). The 
color hypcrfine interaction then leads to attraction in the spin singlet, scalar di- 
quark channel, while the spin triplet, axial vector diquark is disfavoured. Hence 
diquarks should have positive parity and belong to the color 3 representation, and 
so have many properties similar to an antiquark. In lattice QCD this should lead 
to the formation of quark - diquark flux tubes with similar physical characteristics 
to those of quark - antiquark flux tubes. In particular we would expect the long 
range linear part of the quark - diquark potential to have the same slope as that of 
the quark - antiquark potential, corresponding to the flux tubes having the same 
energy density, and we would expect the flux tubes to have similar transverse size. 

The three-quark configurations we use to approximate a quark - diquark string 
are T-shapes, with the origin at the junction of the T. Two quarks are positioned 
one lattice step in opposite directions from the origin (approximating the diquark) , 
and the third is placed from 1 to 12 lattice steps in an orthogonal direction. For 
this part of our work we have used 300 quenched QCD gauge field configurations 
created with the same action as previously. Two hundred of these configurations 
were at at j3 = 4.60 (as in the previous section) and one hundred at (3 — 4.80, to 
investigate the use of a finer lattice. These values of (3 give lattice spacings a of 
0.123 fm and 0.093 fm respectively. 

In Fig. [3] we show an example of the expulsion of vacuum fluctuations and the 
formation of flux-tubes for our quark - diquark configurations. 

The effective potential is obtained from the Wilson loops in the standard man- 
ner: 
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Fig. 4. Comparison of quark - antiquark and quark - diquark effective potentials for /3 = 4.6 (left) 
and /3 = 4.8 (right). 



We obtain stable plateaus for the potentials as a function of r. The quark - antiquark 
potential has the well-known form 

VQG{r) = V -" + *QQr (4) 
where a is the string tension. The three quark potential is EE] 

Vs Q (r) = lv - 1 -j:^- + a 3Q L(r), (5) 

j<k JK 

where Cf = 4/3 and L(r) is a length linking the quarks. As shown in our earlier 
wor L(r) is given by the minimum length of string that connects the three quarks, 
or the sum of distances from the quarks to the Fermat (or Steiner) point. We expect 
that the two string tensions (Jqq and cj^q are equal. In Fig. 2] we plot the extracted 
effective potentials for the quark - diquark and quark - antiquark flux-tubes at each 
of the values of j3 for our gauge configurations. The plots show that our expectation 
is confirmed at both values of (3. Converting length measurements from lattice units 
to fermi we obtain the quark - diquark string tension a^Q = 0.98 ± 0.01 GeV fm _1 , 
which is in excellent agreement with the quark - antiquark string tension Oqq = 0.97 
GeVfrn" 1 . 

We can gain further insight into the properties of the flux tubes by examining 
their profiles close to the quark. We study the values of the correlators C3Q (y) and 
CQq{y) where (y) = (2/1,2/2,0) is constrained to the plane of the color sources, 
and the origin is at the position of either the antiquark or the join of the T. The 
quark is then at the position (£, 0, 0) where £ varies from 2 to 12 lattice steps. First 
we examine the longitudinal profiles of both quark - diquark and quark - antiquark 
flux-tubes along the line (y) = (x, 0, 0) in Fig. [5] As expected, the vacuum expulsion 
close to the diquark is stronger than in the vicinity of the antiquark. However, near 
the quark the two flux tubes show very similar profiles. Similar results are seen at 
(9 = 4.6. 

Next we examine the transverse profiles along a line orthogonal to the midpoint 
of the flux tube, it. along (£/2, y, 0) for £ even, or along ((£ + l)/2, y, 0) for £ odd. 
In Fig. [6] we show profiles of both quark - diquark and quark - antiquark flux-tubes 
for £ = 12. We find that as long as £ is larger than 4 lattice steps, the transverse 
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Fig. 5. Comparison of longitudinal flux tube profiles for quark - diquark (left) and quark - antiquark 
(right) flux tubes at /3 = 4.6 at all longitudinal separations. 




Fig. 6. Transverse profiles for quark - diquark (red lines) and quark - antiquark (black lines) flux 
tubes at (3 = 4.6 (solid lines) and /3 = 4.8 (dashed lines). 



profiles are close to identical. In earlier work ^ we saw that the profiles saturate 
for long enough Euclidean time evolution. Using a fit to a Gaussian profile we find 
that the transverse profiles of quark - diquark and quark - antiquark flux-tubes are 
statistically identical^. 



4. Conclusions 

We have directly compared gluon flux-tubes for quark plus antiquark and three 
quark systems. In the three quark systems we kept two quarks close together (two 
lattice units separation), so that the system approximates a quark - diquark string. 
We found that the string tension in the quark - diquark string was the same as 
for the quark - antiquark string. In addition we compared the vacuum expulsion 
in both sets of flux-tubes. We found that, in the vicinity of the quark, there was 
no measurable difference between the transverse profiles of the quark - diquark 
flux-tubes and the quark - antiquark flux-tubes. 

These findings confirm the expectation from QCD that a diquark has many 
properties in common with an antiquark. In particular the long range color inter- 
action between a diquark and a quark is seen to be the same as that between an 
antiquark and a quark. It would be interesting to repeat this work with dynamical 
quarks, where the effect of screening by quark - antiquark pairs could be observed. 
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